INFRARED BEHAVIOUR OF GLUON AND GHOST PROPAGATORS IN QCD IN 

THE LANDAU GAUGE 
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Non-perturbative formalism of the generalized effective action is used for deriving the Schwinger- 
Dyson equations. In order to clear the domain of integration in the functional integral from gauge 
copies a restriction to the Gribov horizon due to Zwanziger is implemented. In this approach an 
asymptotic behaviour of gluon propagator and propagator of the Faddeev-Popov ghosts at small 
momenta is studied. Such a behaviour is obtained as a result of solving coupled Schwinger-Dyson 
equations in zeroth and first-order approximation. The qualitative agreement of these results with 
the ones obtained before is demonstrated and quantitative difference in some coefficients is found. 



I. INTRODUCTION 

It is widely believed that Quantum Chromodynamics 
(QCD) is the theory which describes strong interactions. 
The physical phenomena at large momenta transfers are 
very well described by the perturbation theory as the cou- 
pling becomes small. Asymptotic freedom allows high en- 
ergy quarks and gluons to be treated as weakly interacting 
particles. This picture, however, starts to break down at 
intermediate momenta and is surely inadequate at ener- 
gies below a few hundred MeV. At such scales the interac- 
tion is strong enough to invalidate the perturbation the- 
ory and one has to employ completely different methods 
to study non-perturbative phenomena like confinement or 
chiral symmetry breaking. One such method is the study 
of the Schwinger-Dyson equations (SDEs). 

Being the non-perturbative equations of motion for 
Green's functions, SDEs plays an important role in this 
sector of theoretical physics Greens's functions which 
can be extracted from these equations, involve an inter- 
esting physics. For example, the non-perturbative be- 
haviour of the infrared gluon and ghost propagators is 
related to the confinement mechanism like the Gribov- 
Zwanziger scenario 2, 3, 15]. Furthermore, certain conse- 
quences related to color confinement can be obtained from 
the infrared behaviour of the gluon propagator in terms of 
unobservable particles, the so-called positivity violation. 

Many attempts have been made to understand the 
gluon propagator behaviour through SDEs. In the late 
seventies Mandelstam initiated the study of the gluon 
SDE in the Landau gauge Neglecting the ghost fields 
contribution and imposing cancellations of certain terms 
in the gluon polarization tensor, he found a highly sin- 
gular gluon propagator in the infrared. This enhanced 
gluon propagator was appraised for many years in the 
literature, firstly because it provided a simple picture of 
quark confinement, since it is possible to derive from it an 
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inter-quark potential that rises linearly with the separa- 
tion, and secondly because the gluon propagator, which is 
singular as l/g 4 , has enough strength to support dynam- 
ical chiral symmetry breaking. However, these results are 
discarded by simulations of QCD on the lattice [(J where 
it is shown that the gluon propagator is probably infrared 
finite. 

Later, infrared finite solutions were also found in 
the Schwinger-Dyson approach [f|. Considering coupled 
gluon and ghost SDEs it was shown that gluon propagator 
is suppressed and ghost propagator is enhanced in the in- 
frared region. Brief overview of the results on the infrared 
properties of the gluon and ghost propagators from SDEs 
existing in today's literature is made in Section HT1 

In the main part of this work (Section 3) the formal- 
ism of the generalized effective action proposed by Corn- 
wall, Jackiw and Tomboulis (CJT) Q is investigated. The 
effective action, which depends on the gluon and ghost 
propagators and the Gribov mass parameter is obtained. 
In Section 4 we construct SDEs for the gluon and ghost 
propagators as well as for the Gribov parameter in one- 
loop approximation. Notice, that the SDEs obtained in 
our work are not exactly the same ones, discussed in lit- 
erature [1, S US EH • Let us clarify this point. 

The starting point in most papers is a full system of 
SDEs. It is a coupled, nonlinear system of equations, 
which contains dressed propagators as well as dressed ver- 
tex functions on the right-hand side. In order to obtain 
a closed system of equations it is necessary to specify 
suitable approximations for these equations. There are 
many truncation schemes, developed for this purpose, for 
example replacing dressed vertices by their tree-level ex- 
pressions. Furthermore, in this formulation, Gribov's pre- 
scription of cutting off the functional integral at the Gri- 
bov horizon does not change the Schwinger-Dyson equa- 
tions, but rather resolves an ambiguity in the solution of 
these equations. 

On the other hand, in our presentation we restrict our- 
selves to two-loop approximation of the effective action 
and obtain a closed system of SDEs without any other 
truncation schemes. These three equations are derived as 
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a result of equating to zero the functional derivatives of 
the generalized effective action with respect to the gluon 
and ghost propagators and the Gribov parameter. Re- 
striction to the (first) Gribov region appears explicitly 
via modification of the free gluon propagator, which en- 
ters the gluon SDE explicitly and, therefore, two other 
SDEs are affected via the dressed gluon propagator too. 

In Section 5 the zeroth-order analysis of the derived 
SDEs is carried out. It is noteworthy that at this order 
we get the results for the gluon propagator and horizon 
condition similar to those previously obtained by Gribov. 
It is shown that the ghost propagator obtained at this 
order and Gribov's one indicate the same qualitative be- 
haviour 1/g 4 in the infrared. However they differ in some 
coefficient. 

The full analysis of the derived one-loop SDEs is de- 
veloped in Section [VII The infrared critical exponents are 
obtained as well as the numerical coefficients for the prop- 
agators. In the last Section (7) the discussion of obtained 
results is presented. 



II. NOWADAYS STATE OF AFFAIRS 

The infrared behaviour of the Landau gauge gluon and 
ghost propagators is an interesting and hot subject. Up to 
now there are two competitive viewpoints on this problem. 

Let us represent the gluon and ghost propagators, re- 
spectively, in the form: 



G%,(q)=C G 



G(q 2 



%1u 



tab 
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where Cg and Cx> are dimensionless constants. We will 
seek for the solutions of the gluon and ghost dressing func- 
tions in the form of simple power laws 



and 



G(q 2 ) OC ( 9 2 ) 



2\ "OC 



(3) 



(4) 



respectively. Coefficients ag and a-p are called infrared 
critical exponents or anomalous dimensions. 

The first point of view is as follows. Using the gluon 
propagator SDE as well as the ghost propagator SDE it 
is claimed that a G + 2«x> = — (4 — d)/2, [&[{j], i.e. for the 
dimension of space-time d — 4 one gets olq + 2a-r> = 0. 
The interesting point here is that self-consistency forces 
an interrelation of the exponents such that they depend 
on one parameter k = a-p = —aojl only. The value of 



the exponent k is in the range 0.5 < k < 1, [§, Itj [ToL ITiT] . 
depending on details of the truncation of the set of SDEs. 
Some of the authors report on the ambiguous result 

aa= —2, olv = 1 and —1.1906, otx>~ 0.5953. 

(5) 

It is worth noticing, that the Kugo-Ojima confinement 
criterion [12[ in terms of infrared exponents requires k > 

0, which means that the ghost propagator should be more 
singular and the gluon less singular than a simple pole. On 
the other hand, Gribov-Zwanziger confinement scenario 
@> 0) E3, EH requires the same condition for the ghost 
propagator, but it needs k > 0.5 for the gluon dressing 
function. Thus, available results for the exponents are in 
good agreement with these two widespread confinement 
scenarios. 

However there is another point of view on this problem. 
Some of the authors [l6[ report that ac + 2ax> ^ and 
correct value of a-p is close to zero while ag is close to 

1. Thus they claim that in the infrared limit the ghost 
dressing function is finite and differs from zero. 

The recent numerical results, i.e. lattice simulations, 
also are not unambiguous [171 ] . 

The aim of the present paper is independently deter- 
mine the infrared critical exponents. 



III. CORNWALL-JACKIW-TOMBOULIS 
FORMALISM 

We start to consider this problem with making use of 
the formalism of the generalized effective action proposed 
by J. Cornwall, R. Jackiw and E. Tomboulis [7]. This 
generalization of the effective action, T(<fi, G), depends 
not only on <fi - a possible expectation value of the quan- 
tum field - but also on G(x, y) - a possible expecta- 
tion value of a bilocal operator T$(x)Q(y) - time ordered 
product of two fields. Physical solutions require 



6<p(x) 



= 0, 



8G(x,y) 



= 0. 



(6) 



In our derivation we identify with the gluon field, 

so G(x,y) is the gluon propagator. We introduce two 
more quantities - the ghost propagator T>(x, y) and the 
so-called Gribov parameter 7. Parameter 7, known also 
as the Gribov mass, characterizes the restriction of the 
domain of integration in the functional integral to the so- 
called Gribov horizon. This restriction is necessary due 
to the existence of the Gribov copies, which imply that 
the Landau condition, d^A^ — 0, does not uniquely fix 
the gauge and equivalent gauge copies still exist in the 
domain of integration of the functional integral 0, [l^] . 

We start from the generating functional for Green's 
functions of nonlocal, composite fields: 



Z(J,K) 



W(J,K) 



(7) 
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Seff =-\[& {F^(x)) 2 + [d 4 xd^C\x)J ab (x iy )C b (y) + 



+ \ [d 4 xd\A;(x)Kf v (x )y )Al(y) + S,(*(x)), (8) 



where 



S z ($(x)) = ld 4 x{ 91 2 r bc (A a ^ b ; + Alvl c ) - A a M)- 



-C d^D^Cf -% c d v {D v ^) ac + uj^d l ,(D v u; ll ) ac + 



+ g(d u cO* c )f abm (D„C)»ip™ + 4 7 4 (A^ - 1)} (9) 

is the expression, which characterizes Zwanziger's formu- 
lation of the Gribov horizon Q (with g being the strong 
coupling constant). It is BRST-invariant together with 
the first term of ©. Functional differential D$ labels the 
product of differentials of all fields entering the integrand 
(effective action (JS])): 



L>$ = dAdCdCdXdipdLpduiduj. 



(10) 



Let us clarify the notations above. Here A is the gluon 
field with components A®. Fields C and C are the 
Faddeev-Popov ghosts which have components C a and 
C . Multiplier A° enforces the constraint d^A^ = 0, char- 
acteristic of the Landau gauge. Fields ip = — ^(tpi+i<p 2 ) 
and 7p = -4= (ipi — itp%) are the pair of boson complex fields 

with components ipf = (p^f = f ahc f b Ji and f abc are the 
structure constants of the SU (N) gauge group. Here we 
use the single index i = ((a, c) for the pair of mute in- 
dices, and i takes on / = d(N 2 — 1) values (d = 4 stands 
for the dimension of space). The fields to and ZJ are the 
Grassmann fields, which have the same components as ip 
and Tp. The vector indices fi, v take on d values, while the 
color indices a, b refer to the adjoint representation of the 
SU (N) group and take on N 2 ~ 1 values. 

We introduced two bilocal sources K^ b (x,y) and 

J ab {x, y) for the gluon field A^ and for the ghosts (7°, C a 
respectively. 

Let G'^' u (x,y) be a possible expectation value of 
TAl(x)A b J (y) and T> ab (x,y) stands for a possible expec- 
tation value of TC a (x)C b (y) : 



G^(x,y) = (TA-{x)Al(y)), 



(11) 



5W(J, K) 
5J ab (x,y) 



= -hV ba (y,x). 



(14) 



The effective action T(G, T>) is defined as a Legendre 
transform of W(J,K): 



Y{G,V) = W{J,K) - 'i ld 4 xd\G ab (x,y)K ab (x,y) + 



+ h d A xd 4 yV ba {y,x)J ab {x,y). 



(15) 



We eliminate K and J in favor of G and T> using the next 
obvious relations: 



2 8T(G,V) 



hSG ab (x,y) 



and 



J ab (x,y) = - 



1 ST(G,V) 



H5V ba (y,x) 
Then according to (|15p. we have 



exp^r(G,P) =exp~W(J,K)> 
h n 



x exp{-^ J d A xd%G;l(x,y)K; b L ,(x,y)+ 



+ i d i xdS J V ba (y,x)J ab (x,y)} 



(16) 



(17) 



(18) 



In order to shorten our notations we will omit dependence 
on coordinates and the corresponding integration where it 
is obvious, keeping it in mind. Taking into account ([7)l. 
© together with (JTHJ) and ([T7J we have 

exp-r (G ,P) =/^e X p-{- (i- ) - K K W AI + 



5G' 



8V ba 



After straightforward computations, for the first exponent 
term in the last expression one gets 



1 



1 



-7 = MiD^XuK - 9f a " c ((dMAlAl) - 



V ab (x,y) 



TC a (x)C b (y) 



From eqs. ([7]) and ((SJ it follows that 

SW(J,K) = l hGab(x } 



(12) 



(13) 
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g * r bc f a P r A b A c A P A r_ 



Here Dp is a free gluon propagator: 



(20) 



(21) 
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The series expansion for r(G, T>) is 

r(G,P)=^rrln(G- 1 )^+^Tr[(^ 1 )^G^]+r 2 (G)- 



- ihTr \n(V' 1 ) ab - ihTr[(S F 1 ) ab V ba ] + T 2 {V) + Const, 

(22) 

where Sf is a free ghost propagator and Dp is a free gluon 
propagator, modified with respect to the Zwanziger term: 



Dp 1 = Dp 1 + 2g 2 1 i N{id 2 ) 



2\-l 



(23) 



The quantity r 2 is given by all two-particle irreducible 
vacuum graphs and it is of order h 2 , because the number 
of loops corresponds to powers of h. The quantities r 2 (G) 
and r 2 (2?) are determined by the expansion by powers of 



(24) 



(25) 



r 2 (G) ^n 2 r ( 2 1 >(G) + o(h 3 ), 



T 2 {V) = h 2 T { 2 1] {V) + (){{?). 

From the expressions (|19p , (j2"0)) the series expansions 
(|2"4"]l and (f2"5|) . using the notation 



T(G,V)=-ihTr \n(V- x ) ab +h?T\ x> {V)-h 2 V 



Here V is the space-time volume of the system and we 
denote the sum of three terms of similar form by means 
of one term: 



Qid 2 Qi = [-C a d 2 C a ~Tpl c d 2 cp« c + uj c l c d 2 ujl c 



(28) 



Thus we have the expression (|2"T)l for the effective action, 
and we have to evaluate the functional integral. But this 
expression includes not only terms, quadratic on the fields, 
which we are able to integrate, but also cubic ones and 
terms of fourth order on the fields. In order to compute 
it, one has to expand cubic and higher-order terms by 
powers of some small parameter in usual way. Following 
Cornwall at al., we will use h (to be precise, we use not H, 
but i/h) as such a small parameter, and this procedure is 
known as loop-expansion. We will retain only terms up to 

h = (^/foj > this corresponds to two-loop approximation. 

In order to select terms by powers h one has to change 
the scale of all 8 fields: 



<E> (x) — > Vft* (x). 



(29) 



After rescaling the fields one can represent the series ex- 
pansion by powers of small parameter h for (|27p in such 
a way: 



SV ba 



T (G,X>) = -in In / Dm^e 



liHTrHG^r + h 2 T^(G) - H 2 GfJ-^?l, (26) 



one obtains apart from the constant ihTrl: 



T'(G,V) = -ih\n /D$exp-{47 4 (iV 2 -l)y+ 

h 



+ -Al{G-'-2g 2 ^N{id 2 )-XvK+' C V 



2 rabcl \a ,J)< 



bc\ 



xfi+VhS^+hs^h^fyfrS^+ns^+oih^J , (30) 

where S q denotes quadratic terms on the fields: 
^-l^(^2.gV^(^ 2 )- 1 )^4-Gl;irL5 F 1 ) Qb G b + 

+ i37 2 / a6c K^ c +^ c )- iAa M)+0«^; ( 31 ) 

The term denotes the cubic ones: 

= i{-gr bc {d^Al)A b ^Ai - gr cb C a d^C b + 



+Qid 2 e^ gf^Xd^DA^l - gr b C\A<fi b 



fact — a o ac b racb — a a \c Jo 

&i o^A^i - gf <Pi dpApCPi - 



„ j'acb — ao ac b racb — ao a c b , / a — ac\ rabra ., 

9 J ^ l duA^ l -g] VidvAtfi+glp&n)} x 



cb—a 



bra . 



x(^)V™ c - -gY'T^AlAtA^Al- hA 



jr^\G) 



SV ba 



" 8G<* 

C b + g 2fab m ^-a^ A X c \ blf rnc_ 



A 



(27) 



+ g{d v U5l c )f abm {d v C) b vTl 



(32) 



and 



= i{--g 2 r hc r r A\AlAlAl - A 



q gMG) b 

" 8G± v 



l^ST * + 9 2 f abm (d^)AtC Xb ^ c ) (33) 



C 
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denotes the terms of the fourth order on the fields. 

Carrying the term with 7 in front of the sign of the 
integral in eq. ([30]) , then collecting terms by powers of H 
in the last multiplier of this equation and retaining only 
terms of order H, we obtain the two-loop contribution to 
the effective action of the theory : 

T'{G,V) = 4 7 4 (7V 2 - \)V- 



=\9 2 r bc ) 2 {-G, P (x, y)(d^G va {x, y))(d^G^(x, y)) + 

+ G pp (x,y) (d^G va {x,y)) (d^G v<7 {x,y))}. (38) 

For four gluon fields at the same space-time point A^(x), 
one gets: 



-ihln J L><M 4 e s 
Expanding in h we obtain 



(34) 



T (G,£>) -4 7 4 (A^ -1)V = 



= -ihln / L><M 4 e s 



~ih 2 (sW + 1(S(^) 2 \, (35) 



The first term is the result of integration e s ' on all 8 
fields (see the notation (fTU)0 . together with the notations 
(|3ip and (|28[) . It gives, omitting unimportant constants 
— iMnh 4 and \KTr In [d^d^]: 



ihln J D^h 4 e s " = -ihln [Det(V- v ) ab ] 



(36) 



Note that integration over the fields (p,Tp and u>,TD gives 
unity as well as over A and A. 

The second term of (|35|) is a sum of all vacuum ex- 
pectation values of all the field configurations, included 
inside brackets ( ) . One has to write down all these aver- 
ages explicitly, making possible use of the Wick theorem 
for the computing every vacuum expectation value. Thus, 
all configurations with non-zero expectation values, which 
contribute to the second term of eq. ([33)) are written as 
follows: 

S m + l(S^) 2 ^ = - l -g 2 f abc f apr (AlAlAVAl) - 



g 2 f a bcfP rs (( d ^)A b At(d p Al)A r p Al) 



- -/r ch f rs (e a (d p A^e b )e p (d„A:e s )) . (37) 

Let us consider these expectation values, making use of 
the Wick theorem. Here we write down only the result of 
the computations. 

1 ^/^^A^ll^)^)^) 



~/f abc f apr (Al(x)AUx)AP(x)Ai(x)) = i .g 2 (r hc ) 2 x 



x {G^(x, x)G flv (x, x) - G^(x, x)G vu {x, x)}. (39) 

For anti-commuting ghosts (the first three fields are lo- 
cated in the same space-time point, say x, and let y be 
space-time point for the location of the last three fields), 
one gets: 



-\g 2 r b F rs ({d\r ] c a )A^c\ x ^c p )Alc s 



(v) 



= -\ 9 2 {r b ) 2 {d^V{y^)d^V{x,y)G, u {x,y)}. (40) 

Now with a close analogy to the ghosts, for Grassmann 
fields lu,ZJ we get 



~9 2 f acb f prs ( (d^uf t )Alu\\ [x) {dWufAAW | 



= ^9 2 {f acb ) 2 {d^G%{y, x)d^G% (x, y)G fiU (x, y)}. 

(41) 

Here Gfj(x, y) is a tree level propagator for the fields to, uj 
and indices i,j take on d(N 2 — 1) values. For a pair of 
complex conjugate bosonic fields (p,Tp we get 

1 2 f a c bfP rs / ( ^)-a } ^ ( ^ fyy)^ A> s ] 



= -l9 2 (f acb ) 2 {di x) Gf J (x,y)di^Gf J (x,y)G^,y)}, 

(42) 

where Gf^ (x, y) is a tree level propagator for the fields <p, ip 
and indices i, j are the same ones as in the previous case. 

Finally, for the quantity T (G, V) from expressions (|35|) 
- (O one has: 



r'(G,X>) = -^Trl + 4 7 4 (7V 2 - l)V - ihTrln{V~ 1 ) ab - 



klsnab jTj \G) j_ T,ba jTg JP) , i <2( rabc\2^ 
n V^IXU snnh h U 5V ba b 19 (J ) X 
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x {G p „G pi/ - G^Guu) + \g 2 (f acb ) 2 [~d^Vd^VG^- 



lh lw¥ ( TrlnP_1 ^ +TrS F 1 ( P )V(p) - 1) 



x [-G„ p {d p G„ a )(d p G pr7 ) + Gp P {d p G va ){d p G va )]}. (43) 

As it follows from the explicit definition of the propagators 
Gfj and Gfj (it can be found during the explicit integra- 
tion of the expression (|36|) for S q with the notations (|3"Tj) 
and ([25])). = G?-. Inasmuch as these propagators en- 
ter the eq. (I43|) via terms of equivalent form, their whole 
contribution cancels: 

- d^G%d^G%G pv + d^G%d^G%G pv = 0. (44) 

Thus, the effective action has no dependence on auxiliary 
unphysical fields tp, Tp and ui, ID. 

Finally, taking into account that for the gauge group 
SU{N): 



rade rbde Jtf§ a ^ 



(45) 



the explicit expression (22) for the effective action of the 
theory reads as: 



T(G,V) = 4 7 4 (V 2 -1)V + -ihTr[ HG~ l )±+ 



-{D-^G^-^-ihTr [ln(p- 1 ) afc +(5 F 1 ) Qb I? ba -l] - 

■^ih 2 g 2 N5 cc {^{G^G^-G^G vv )-d^Vd^VG lxv - 
- G up (d l j l G V(T )(dpG ^ + G pp (d p G^)(d p G„ a )}. (46) 



It is convenient to rewrite the expression (J46J) in momen- 
tum space, using the Fourier-transformed propagators, 
defined as follows: 



G(p) 



d 4 xe^ p{x - y) G(x~y); 



(47) 



(Inasmuch as we seek for the ground state, we consider 
the case of translation-invariant solutions, i.e., we take 
G(x, y) to be a function only of x — y, and similarly for 
the other propagators.) 

Passing to momentum space, for the eq. (|46|) . divided 
by the volume V, we get 

~f(G,X>) = 4 7 4 (V 2 -1) + 



'Ml^W (TrlnG-\p)+TrD^( P )G(p) - 1 



(27rfi) 



y 1 {2irh)» 



%{G P v{p)G pv (r)-G pp {p)G vlJ {r)}~ 



{V{p)V{r)G pu (p-r)p p r v + 



+G vp (p)G va (p + r)G p<y (r)(p + r) p r p - 



- G pp (p)G va (p + r)G ua (r)(p + r)^}. (48) 

Thus we have the two-loop generalized effective action, 
written in the momentum representation. 



IV. SCHWINGER-DYSON EQUATIONS 



In order to obtain a closed system of equations on the 
gluon propagator G^ v {p), the ghost propagator V ab (p) 
and the Gribov parameter 7 one has to take variational 
derivatives of the effective action (|4*51) with respect to 
these parameters and equate each of them to zero (see 
eq. ©)• 

Before doing this, it is worth to remember that in the 
momentum representation the expression (|23p for the free 
boson propagator modified with respect to the Zwanziger 
term is given by 



D F \p) = D F \p) + 2g^N%. 

p z 

From here on we will put H = 1. 

Varying expression (|4*5|) with respect to 7, we get 

1 5f(G,V) 
V <?7 



(49) 



4 7 d 



A(N 2 - 1) ~g 2 N 



d 4 P m G»£(p) 



(27T) 



rTr 



P 



(50) 



and after rejecting the trivial solution 7 = 0, it looks as 

Ng 2 f d^p „_G%lp) 



1 



A{N 2 -I) J (2tt) 



rTr- 



V 



(51) 



Variation with respect to T> ab (q) is given by 
1 ST(G,V) 



= o = i(T>- 1 {q)-Sp 1 (q)) ab + 



V SV ab (q) 



7 



-g 2 N8 



ab 



d 4 r 



{G-\q))^={b- F %)) a l-9 2 N5 a 



d 4 p 



T>(p)Gfj, v (p- q)p fl q v ] 



(52) 



Due to the property G(r) = G(—r) this expression after 
changing of mute indices p — > r and /i — > v in the last 
term, can be written in the form 

(V-\q)) ab = (S F \ q )) ab + 



2 cab 



Ng 5' 



d 4 r 



(53) 



Variation with respect to (g) is given by 



1 6T(G,T>) 
V SGfM 



= {) = U-G-\ q ) + D F \ q ) 



ab 



WNS ab x 



V(p)V(p - q)p fl {p - q) u + (gluon loops). (54) 



Some simplifications can be done in the expressions above. 
In the eq. (|5"Tj) we take the trace on color and Lorentz 
indices. In the Landau gauge the gluon propagator is 
transverse on vector indices: 



Gf v (p) = G(p 2 )P^(p)6 ab 



(55) 



where 



d 4 p 



{D{p + q)V(p)(p + q) ll p v } + (gluon loops). (59) 



It is worth noticing, that the Gribov parameter 7 does not 
enter its equation of motion (|57| explicitly, but through 
the modified free gluon propagator Dp(q) (according to 
(|49p ). which is present in equation (|59[) for the gluon prop- 
agator. 

We wish to determine the asymptotic form of the prop- 
agators at low momentum. For this purpose we let the 
external momentum in the SDEs written above be asymp- 
totically small. In this case the loop integration will be 
dominated by asymptotically small loop momentum, so 
the propagators inside the integrals may also be replaced 
by their asymptotic values. As one can see from Section 
2, where it is argued that the ghost propagator should 
be more singular and the gluon less singular than a sim- 
ple pole, gluon loops integration at low momentum in eq. 
(|59"|) may be neglected. 

Let us consider the eq. (|55|) . One of the form of the 
horizon conditions, that guaranties the absence of the Gri- 
bov copies, is written as [l8| 



lim [q 2 V(q 2 )] 1 = 0. 



To impose this condition, we divide eq. (|58p by q 2 , and 
obtain (after the factorization of color indices) 



= i + Ng z 



d 4 



(M 



-V{r 2 )G(r 2 ) [1 - cos(gr)] . (60) 



(56) Subtracting this equation from the previous one l|58[). we 



is the transverse projector. In the eq. (|53")) we shift the in- 
tegration variable according to r — > r+q and use eq. (|55[) . 
The similar shift of integration variable will be applied to 
the eq. ([53]). 

Now we are ready to write down a closed system of 
equations of motion for the Gribov parameter, the ghost 
and gluon propagators: 



1 = iNg J 



d A p G{p) 
(2tt) 4 p 2 ' 



(57) 



V-\q 2 ^ 



-Ng' 



d 4 r 



1 2 

r q 



G(r 2 )[V(r 2 )- 



(qrf 



(61) 



We now turn to the SDE for the gluon propagator (159[) 
with comments below it. We apply the transverse projec- 
tor and take the trace on color and Lorentz indices and 
obtain taking into account eq. (f49|) : 



{V-\q)Y b = {S F \q)) ah + 



G-\q 2 ) = i q 2 + 2Ng 2 1 ^- 



2 cab 



Ng z 5- 



(27T)' 



-V(q + r)G(r 2 ) 



r 2 q 2 — {qry 



■ (58) 



p 2 q 2 



(pqf 



(62) 
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Notice, that in the limit q — > the first term vanishes. 

In order to simplify further analysis and compare our 
results with other authors, we rewrite SDEs for the gluon 
and ghost propagators in Euclidean space: 



V- 1 (q 2 )=Ng 



d 4 p 
(2tt) 4 



G{p 2 )[V(p 2 )-V{{q+p) 2 )} x 



p 2 q 2 — {pq) 2 



1, -a 94 1 Ng 2 fd 4 p 
G- 1 {q 2 ) = 2g 2 1 A N-+ " ' 



Q 



3 J (2tt) 



p 2 q 2 — {pq) 2 



(63) 



V{{p + q) 2 )V{p 2 )x 



(64) 



This system of SDEs for the gluon and ghost propaga- 
tors in Landau gauge, written in Euclidean space, will be 
analyzed in next sections. 



V. ZEROTH-ORDER ANALYSIS OF SDES 

Now we consider the first approximation which corre- 
sponds to retaining only the first term in the right-hand 
side of ([55]) without any changes in other equations. It is 
in fact the zeroth-order approximation. In eq. (|59[) the 
contribution at this order gives 



/ , « \ ab 

G (1) (<Z) = 



2Ng 2 j 



.4 I 9n* 



(65) 



To see this one has to inverse left-hand side and the first 
term in the right-hand side of eq. (|59p and rewrite the 
result in Euclidean space. 

After substitution of this approximation to the Eu- 
clidean form of eq. (|57|) . we get 



;Ng' 



d 4 q 



1 



(2tt) 4 q 4 + 2g 2 Nj 4 



(66) 



Note, that the gluon propagator (f65|) and the gap equation 
(|6"6"|) were first derived by Gribov 2]. Furthermore, so- 
called horizon condition (|66[) was obtained by making use 
of the zeroth-order approximation (|65|) . so we see that 
SDE for the Gribov parameter (I57|) defines full horizon 
condition. 

We now turn to the ghost SDE and apply zeroth-order 
approximation (|65[) to it. For this purpose it is more con- 
venient to use Euclidean form of eq. (|53p . Substituting 
the explicit expression for G llv {q — r) (according to (|65p ) 
in the horizon condition (|66p one gets 



V- 1 {q)=Ng 2 q ti q v 



1 



(2tt) 4 r 4 + A 4 



[l - r 2 V(q - r)] 



(67) 



with A 4 = 2Ng 2 1 4 . 

This is a nonlinear integral equation for the T> - func- 
tion and we now have to solve it. One of the possible 
ways to do this is iteration scheme with a free ghost 
propagator Sp{q) — I/9 2 as a zeroth-order approxima- 
tion 2?(o)(<z) °f 2? (<?)■ I n other words, in order to obtain 
first-order approximation D^{q) of r D^ 1 {q), we have to 

replace T>(q — r) — * D(o)(<7 — r) = l/(q — r) 2 in the right- 
hand side of eq. (|6T|) . Doing this we get 



TT^q) =Ng\q v 



d 4 r 



1 



(27T) 



-A 4 



q 2 ~2qr 

¥- 



(68) 

Notice, that exactly the same expression was used by 
Gribov in order to obtain the infrared behaviour of the 
ghost propagator 0. 

We emphasize that while Gribov in fact solved the equa- 
tion for the ghost propagator approximately, we derive 
some integral equation and we look for the exact solutions 
to it. 

After integration of the right-hand side of (|68| for the 
first-order approximation of T>(q) one gets [l9j: 



%)(<?) = 



128ttA 2 1 
3g 2 N q 4 



(69) 



Notice, that further making use of such iterational scheme 
for the eq. (f67|) results in incompatible system of equa- 
tions. Thus the natural question appears: is the itera- 
tional procedure, used for solving the equation (|6T|) . ad- 
missible? The obvious answer is " no" . Because this itera- 
tion procedure is not converging. Consequently, perturba- 
tion theory does not work and the result, obtained above, 
is in general incorrect, though it gives qualitatively right 
behaviour. In order to solve the equation (|6"T)) correctly 
some another method is needed. 

Let us try to solve the equation for the ghost propagator 
(|67p , written in the form, similar to (|63p (with the explicit 
gluon propagator (|65|0 in the infrared: 

[gV{q 2 )] 1 =N 9 J^- [V{r 2 ) -V{{q + r) 2 )] x 



\r 2 q 2 



{qrf 



A 4 



(70) 



In order to solve this equation, one can seek for the solu- 
tion in the form 



9V{q 2 



A 



(<7 2 ) 



l+a ■ 



(71) 
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Here A is some coefficient and a is called infrared critical 
exponent. After evaluation of the integrals, one gets 



2 \4 



A = 



32?HA 
N 



(72) 



Thus, the exact solution of the equation ([67]) . which de- 
scribes the infrared behaviour of the ghost propagator is 



V(q 2 ) = 4^^1l. 

gVN<l 4 



(73) 



Notice, that the previously obtained solution ([69)) quali- 
tatively coincides with ours, but has another coefficient. 
Thus our zeroth-order analysis qualitatively confirmed 
that in the infrared region the gluon propagator is sup- 
pressed (|65|) . while the ghost propagator is enhanced ([73)) . 



VI. FULL ANALYSIS OF THE DERIVED 
ONE-LOOP SDES 



We now turn to the SDEs ([53"]) and ([M]) . One can solve 
these equations using "weak-angle-dependence" approxi- 
mation for the ghosts: 

V((q + p) 2 ) ^V(q 2 )Q(q 2 -p 2 )+V(p 2 )Q(p 2 -q 2 ). (74) 

Here 0(x) is the ordinary step function. Substituting this 
approximate expression to the ghost SDE (|63p and gluon 
SDE ([Ml) and then carrying out the angles integration, 
one gets 

rl 2 

V-\q 2 ) = Aq 2 / dp 2 p 2 G(p 2 ) [D(p 2 ) - V(q 2 )] , (75) 
Jo 



G-\q 2 )= ^+B fdp 2 p i V{p 2 )V(q 2 )+[dp 2 p i T> 2 ^) 

q Jo J„2 



with 



A 



_ 3Ng 2 



B 



N9 u A^2iV,V. 



16(2tt) 2 ' 16(2tt) 2 ' 
Solutions may be found in the form of power-laws: 
G(q 2 ) = C G (q 2 )- a , 



V(q 2 ) = Cvtf) 



2\-/3 



(76) 
(77) 

(78) 
(79) 



with Cg and Ct> - dimensionless constants, a and (3 - 
some exponents. It turned out, that straightforward sub- 
stitution of anzatz 1(78]). (|79| in eqs. ([75")). ((76)) suffers 
from certain mathematical difficulties and it is more con- 
venient to do first some transformations in the equations 



above. We divide equation (fTS"|) by q 2 then carry out dif- 
ferentiation with respect to q 2 in derived equation and in 
eq. (f76|) . and obtain: 



2 

\v-Hq 2 ) + ^V- 2 {q 2 )v\q 2 )=AV\q 2 )fdp 2 p 2 G{p 2 

q 4 q 2 Jo 



(80) 



G- 2 {q 2 )G\q 2 ) = ^-BV(q<)ldp'p*V(p 2 ). (81) 



A 4 



„2 „4 T1 /_2\ 



Now the anzatz (|78j) . ([79]) can be applied. We substitute 
this power-laws into our equations and carry out all the 
integrations (under the assumption of /3 ^ 3 and a^2). 
Then we get 



H3 + i)^vr 2 + ^J 



(q 2 )- ~ a+1 = 0; (82) 



-aCa\q 2 ) 



A 4(^-2 



(q 2 )- 2f}+2 = 0. (83) 



We equate exponents of q 2 in both sides of this equa- 
tions and obtain following equations for the parameters a 
and (3: 

P-2 = -P-a + l (84) 



a- l = -2 = -2/3 + 2, 
which give us a unique solution 

a = -l (3 = 2. 
In terms of the infrared critical exponents 
aa = ot — 1, cud =0—1, 
(see ©, g])), one has 

Qg = —2, dp = 1. 



(85) 
(86) 
(87) 



We now turn to the equations ([82)1 and ((83)) in order to 
obtain the coefficients Cq and C-p. Thus we get 



a 



t> t -AC-pCc — 0; 



A 4 - 2BC|, = 0. 



Along with the notations ([77)) this gives 

Cg = ±[1 + 3BA- 1 ]=1L; 



Cr 



3 -A-^ = A^ 
2 G 



gVN 



(89) 



(90) 



(91) 



(92) 
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with A 4 = 2Ng 2 j 4 , where 7 is the Gribov parameter. 

Consequently, the solution of the SDEs ((55)) and 
is the following infrared behaviour of the gluon and ghost 
propagators, respectively: 



G(1 2 ) = 



V(q 2 ) = 4 — =—r. 



(93) 



(94) 



Thus one can see that the onedoop results are in good 
qualitative agreement with the tree-level ones. The coef- 
ficients (|91|) . (|92|) obtained at this order, differ from the 
zeroth-order ones (|65|). (|73j) by numerical multiplier and, 
hopefully, are more precise. 

As can be checked by numerical computation, the ap- 
proximation (|74[) have no influence on the qualitative be- 
haviour of the propagators, but rather gives some addi- 
tional constant terms, unimportant for the asymptotic. 



VII. DISCUSSION AND CONCLUSIONS 

Formalism of the generalized effective action with im- 
plemented restriction to the Gribov horizon is used to 
investigate the infrared behaviour of the gluon and ghost 
propagators. The system of the SDEs on the gluon prop- 
agator ([59]) . the ghost propagator ([58]) and the Gribov 
parameter (|57|) is derived in this approach. These equa- 
tions are obtained using two-loop approximation for the 
generalized effective action. 

The SDE for the Gribov parameter is in fact the full 
horizon condition, which reduces to the Gribov's one at 
zeroth order. Furthermore, the gluon propagator in this 
(zeroth-order) approximation is exactly the same expres- 
sion, first pointed out by Gribov. The expression for the 



ghost propagator in zeroth-order approximation agrees 
qualitatively with the Gribov's one but has different nu- 
merical coefficient. 

The obtained SDEs are also solved in the first-order 
approximation. The infrared critical exponents and coef- 
ficients for the gluon and ghost propagators are obtained. 
Thus, we derive the expressions, which describe the in- 
frared behaviour of the gluon (|9"3"|) and ghost (|9~4"|) propaga- 
tors. The solution for the infrared exponents was obtained 
before whereas the coefficients for the propagators at 
this order are derived first in the present work. 

As expected, the first-order results agree with the 
zeroth-order ones with different numerical coefficients at 
small momenta. For the gluon propagator this difference 
reduces to 1/2 and for the ghost one - to y/2. 

Thus, we can see that cutting off the measure of the 
functional integral at the (first) Gribov region deeply 
modifies the structure of the propagators. This leads to 
infrared enhancement of the ghost propagator and to sup- 
pression of the gluon one, as can be seen from our calcu- 
lations, described above. It is significant, that such be- 
haviour of the pro pag ators supports the picture of color 
confinement, [2t l3l. Il5j. which we believe to occur. 

It would be interest ing to study the possibility of gener- 
ating the gluon mass [22( in the effective action formalism 
of CJT when the Zwanziger horizon condition is imple- 
mented [20j . 
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